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Abstract 


A formalism is developed for evaluating the momentum distribution 
for proton production in nuclear abrasion during heavy ion collisions 
using the Glauber multiple- scattering series. Several models for the one- 
body density matrix of nuclei are considered for performing numerical 
calculations. Calculations for the momentum distribution of protons 
in abrasion are compared with experimental data for inclusive proton 
production. 

Introduction 


The prediction of the secondary distribution of 
particles produced in heavy ion collisions is re- 
quired for performing cosmic-ray transport calcula- 
tions (ref. 1). The fragmentation cross sections in 
heavy ion reactions may be described in the abrasion- 
ablation models in which abrasion describes the 
knockout of nucleons and clusters during the projec- 
tile and target overlap and ablation describe the de- 
excitation of the projectile or target remnants. The 
abrasion cross sections are typically described us- 
ing a geometric model (ref. 2) or semiclassical meth- 
ods based on the Glauber or eikonal approximation 
(refs. 3-5). 

In this report we consider the use of the Glauber 
multiple-scattering model for formulating the pro- 
ton momentum distribution from abrasion. This for- 
malism will be useful for obtaining both a better 
understanding of the physics of nuclear abrasion and 
estimates of secondary yields of light particles from 
direct processes. We will also relate the inclusive nu- 
cleon scattering observables to the internal nuclear 
density matrix that differs from the intranuclear cas- 
cade (ref. 6) or hydrodynamic model (ref. 7) descrip- 
tions of these observables. The calculations may then 
provide tests on nuclear structure calculations for the 
internal momentum distribution when all secondary 
nucleon mechanisms are considered. The same for- 
malism describes proton or neutron production from 
nuclear abrasion; however, because of increased mul- 
tiplicity for evaporation neutrons as compared with 
that of protons, we present calculations for proton 
production in this report. 

Several mechanisms for proton (or neutron) pro- 
duction occur in heavy ion collisions. Figure 1 il- 
lustrates the mechanisms for proton production in 
nucleon-induced reactions. The contributions from 
the scattering of the incident nucleon from the tar- 
get ground and low-lying excited states are close to 
the beam energy. Below the beam energy we find, 
first, the quasi-elastic peak where the projectile has 
knocked target nucleons into continuum states, and 
this is followed by the nucleon excitation peak which 



Figure 1. Schematic of mechanisms of proton production in 
nucleon-induced reactions. 


at intermediate energies is dominated by the delta 
resonance. The nucleon excitation peak is observed 
to grow in importance with increasing beam energy 
and with the production angle of nucleons. For inci- 
dent neutrons (protons), a quasi-elastic peak of pro- 
tons (neutrons) is observed because of charge ex- 
change. At lower energies, the two contributions to 
proton production are (1) the evaporated particles 
from the deexcitation of the target nucleus and (2) a 
higher energy contribution from the knockout or cas- 
cade particles and the multiple scattering of such par- 
ticles inside the target. The cascade particles extend 
out to several hundred MeV because of the multiple 
scattering and the internal Fermi motion of the tar- 
get, and these particles will overlap with the quasi- 
elastic peak at forward angles if the beam energy is 
not too high. 



For heavy ion reactions, many of the same mech- 
anisms will contribute to proton production however 
modified, because both the projectile and target nu- 
clei are a source of evaporation and knockout pro- 
tons. An illustration of the contributions to proton 
production in heavy ion collisions is shown in fig- 
ure 2. The internal Fermi motion of the projectile 
allows secondary protons (or neutrons) to be pro- 
duced at appreciably higher energies than the beam. 
The nucleon excitation peak will still appear in the 
intermediate region between the projectile and tar- 
get. One important difference between the nucleon- 
induced and heavy-ion-induced proton production is 
the possibility of coherence effects which may lead to 
stronger multiple-scattering contributions for heavy 
ion reactions. In this report we will describe the cal- 
culation of proton spectra from the knockout mecha- 
nism. The contributions from nucleon excitation and 
evaporation will be described elsewhere. 



Figure 2. Schematic of mechanisms of proton production in 
heavy ion reactions. 

Symbols 

A mass number 

b impact parameter, fm 

c speed of light, m/sec 


E energy, MeV 

F fragment 

F* prefragment 

/ scattering amplitude, fm 

/nn nucleon-nucleon scattering 

amplitude 

Im imaginary part of function 

K initial wave number, fm -1 

k wave number of knockouts, fm -1 

M mass, MeV/c 2 

myj nucleon mass, MeV /c 2 

N 1 single-collision term 

n number of abraded nucleons 

n(p) momentum distribution 

P projectile 

p momentum, MeV/c 

Q defined in equation (13) 

q momentum transfer, fm~ ] 

r internal nuclear coordinate vector 

S n separation energy, MeV 

s transverse part of r vector 

T target 

Tp* kinetic energy of prefragment, MeV 

t time, sec 

X final target state 

F profile function 

S Dirac delta function 

A n defined in equation (20) 

£ defined in equation (22) 

p(r, r') one-body density matrix, fm -3 
p{r) one-body density, fm -3 

cr cross section, mb 

X eikonal phase 

D eikonal inelastic collision term 
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Abrasion Theory 

In the Glauber model the scattering operator for nucleus-nucleus collisions is written as 


f=~ J d 2 be^ h V(b) ( 1 ) 

(refs. 8 and 9) where K is the projectile-target relative wave number, b is the impact parameter, q is the 
momentum transfer, and the profile function representing the multiple-scattering series at high energies is 

r(b) - 1 - ni 1 - r ^'( b - 8a - Bj)] ( 2 ) 

a,j 

where a and j label the target and projectile constituents, respectively. In equation (2), F aj is the two-body 
profile function with the internal coordinate having components r — (s ,z). 

The scattering amplitude of equation (1) is related to the production cross section for a projectile nucleon 
from the abrasion process by 


da 

dk 




(TP|r t (b / )|XF :,: k i )(k i F*A|r(b)|PT) 


( 3 ) 


where k j denotes the wave numbers of the abraded nucleons, F* denotes the prefragment (with A F * = Ap-n), 
and we have inserted initial and final states in equation (3). 

The state dependence of the final target energy prevents closure on these states from being automatic in 
equation (3), although when energy conservation is not considered it is made outright (ref. 3). If we consider 
the change in energy of the target from the collision 


e T ~ E x = E t - yj (p T - q) 2 + Mj, 

( 4 - M 

= Et 1 - A 1 + 




q 2 + M l - M 2 


X 

E 2 


( 4 ) 


where M T and M x are the mass of the target in the initial and final states, respectively, we expect that 
performing closure on |X) will be valid for sufficiently large values of E F . In proceeding with closure on the 
final target states, we now write 


where we define 


-j 

£ = (2^/ 


.(b,b / ,k i ,q,£? F *)=(r|| J ft 7^3 6(Ei-E f ) 


x (P|rt(b / )|F*k i )(k J -F*|T(b)|P) >|T) 


( 5 ) 


( 6 ) 


3 



In order to consider the energy-conserving delta function in equation (6), we introduce the Fourier transform 
pair 

Vn{t) = J dEe~ lBt a n (E) ( 7 ) 

and 

ME) = J £ e iEt a n (t) ( 8 ) 

In the projectile rest frame we have 


n i_2 


E{ E j — S n Tp* -j- Ep — Ex — ^ ^ 

j-l 2m iV 


(9) 


where S n is the separation energy and Tp * is the recoil energy of the prefragment, including any excitation 
energy of the prefragments. We next go into temporal space in order to consider the dkj integrals in equation (6) 
and thus rewrite equation (9) as 

k? 


Ei 


~ E f = E~Y. 


3=2 


2m tv 


( 10 ) 


From equations (6) and (7) by using equation (10), we find 


a n (t) = <u / n 


dk. n 


3 = 2 


L(2tt) 3 J 




x (P|rt(b')|F*k J XF*k,-|r(b)|P)l|T) 


(ii) 


In order to simplify equation (11), we first factor the profile function into projectile participant and spectator 
terms as 


where 


r (b) = 1 


Ap n 


n Q/( b — s z ) jj Qj(b - Sj) 

l=n + 1 j — 1 

(12) 

Aj’ 

Qj = II( 1 ~ r «j) 

(13) 


a—X 


In the abrasion model the orbits of the prefragments are assumed to be nearly the same as those of the 
projectile. This is consistent with the use of the impulse and frozen nucleus approximations at high energies. 
A completely factored form in the participant and spectator coordinates is assumed for the projectile wave 
function 


! P ) = \F)\<f>n) (14) 

where |F) and [ <j> n ) are the wave functions of the core (spectators) and of the knockouts (spectations), 
respectively. The antisymmeterization is ignored in equation (14), which should be accurate if the mass of 
F is much larger than the knockouts. Antisymmeterization in the subsystems of |F) and \<p n ) may still be 
included. By using plane-wave states for (ky) and substituting equations (12) and (14) into equation (11), we 
find that 


n(t) = (T\ j ( A ?)(F\ n Qj (W - s') \F*)(F*\ JJ Q z (b - s t )|F) J dr, dr[ q\ (b' - sj) Q( b - Sl ) 


X II 

3 = 2 L 


dki 

(2^ 


; d F dr'j e* k i x t e Q\ (b' - s') Q ;/ (b - Sj ) 


0 n (ri,...,r n )4(ri,...,r;) l]T) (15) 
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where x,- - rj - r'-. Using the coherent approximation for the target wave function in the intermediate states 
and the independent particle model for the fragment wave function leads to 

e n (t) = (^f) P^- n ( b, b') A„_!(b, b', t ) (16) 

where the function P(b, b ; ) describes the projectile spectators as given by 

VAP ~ U ( b. b') = (TF\ n Q; (b' - s|)|F*><F*| <J,(b - s,)| FT) (17) 

l l 

We next perform closure on the prefragment states in equation (17) because we do not consider coincidences 
with individual states. After closure we find that 

P^(b,b') = {TF\ JJ Qj (b' — s;) Qi(h - s;)|FT) (18) 


In equation (16) we have defined 


5k" = (2^ / **' e W'(ry)<2, t (b'-s')<3;(b-s) 


(19) 


where p(r,r' ] is the one-body density matrix of the projectile given by p(r, r') = MrjdMr'). Next, from 
equation (15), after evaluation of the integrals over k, for j > 2, we find 

A n -l(b, b', t) = (T\jf[ (dr J dr'(^) 3 / 2 e- m ' VI f / 2 i V(r J -, 

J 2 

xQj(b' -s'^Q^b-Sj) 


|T) 


In energy space, 

An-! (b, b', Tp * , k) = (Tj J 


3 = 2 


m N ( 1 \3(n-l)/2 ^ 1 ) /2 I- 1 


2 \2n ) 


[3(n-l)/2]-l [3(n— 1)/2] 

x n~ 1 


J \Z(n - 1 )/2] - 1 1 X n~ 1 ) I : T ) 


( 20 ) 


( 21 ) 


where J$ is the cylindrical Bessel function of the first kind of order m and where 


£n- 1 


2m n 


Tp* -f- 


k 2 

2m at 


~S n 


(E r - E_x ) 


(22) 


and 


n 



(23) 
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For n — 1, we have Aq — 6(Tp*). If we assume forward-peaked density matrices (about x — 0), a small 
argument expansion of the Bessel functions can be developed (refs. 10 and 11) that results in 


A n _i(b,b / ,r F *,k)«C , n -.i 

xA?- 1 (b,b',-|^) + 0« 4 x 4 ) 

\ vra - 1/ 


TF‘ + ~-S n -(E T -E x ) 


n — 1 


(24) 


where, for example, C\ = 1, C 2 = tt/4, C 3 = tt/ 105, and C 4 - tt 2 /240. 

The nucleon momentum distribution from abrasion then takes the form 


d<J 

dk ^ 


A P \ 1 


dNi 

~dk 


)(2^j2 / ^didhdh' V Ap ~ n (b,b') 

J d7>. A„_i(b,b',T F .,k) 


If 

we may approximate equation (25) as 


Q ^ + Mp — M.\ <C Ep 


da 

dk 



d 2 bV Ap ~ n { b) — f 
dk J 


dTp* 


A n -i(b, Tp*, k) 


(25) 


(26) 


(27) 


The result of equation (27) suggests that for Ap 1, the momentum dependence of higher production terms 
in > 1) should be similar to the leading-order terms as described in the appendix. This supports the success 
of the hard-scattering model of Hatch and Koonin (ref. 12) where only the single-scattering mechanism was 
used to predict the shape of the inclusive proton distribution in heavy ion collisions. The model developed 
herein differs from the hard-scattering model by our use of the target closure approximation in which the effects 
of smearing of the secondary momentum from the target knockouts are not considered but are replaced by 
averages represented by the target density and by using only on-shell two-body amplitudes. Also, the Glauber 
model has a much fuller multiple-scattering structure than the hard-scattering model. 

In order to include the effects of final-state interactions (FSI) of the nucleon knockouts, we use the eikonal 
model described in reference 11 in which the plane waves are replaced by the distorted waves for the nucleon- 
projectile recoil interaction evaluated at the relative energy between the knockout and recoil. Modifying 
equation (19) as in reference 11 gives 

(^ 3 / *^e‘ k *p(r,rOe- ata x ( - , (»)ol (b '- 8 ' ) <? l Cb-s ) (28) 

where x ^ is the outgoing eikonal phase. Equation (28) ignores off-shell effects, whereas the energy dependence 
of the FSI is included and a medium modified interaction is assumed as described in reference 11. 


Optical Limit for Profile Functions 

For ApAp 1, the optical limit of the profile functions occurring in the previous expressions may be used 
(refs. 8 and 13). From reference 13 we find in the optical limit that 


V Ap (b,W) 


— exp 
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(29) 



where the eikonal phase is 


*( b ) ~ 27dfe N ^ J ' ^ e * q b ^V(0) Ft{q) /nn(?) (30) 

with F denoting the one-body form factor and /nn denoting the two-body scattering amplitude, which we 
represent by 

fm = °MN( Q + < )*NN e -l^ (31) 

where ctnn is the two-body total cross section, B is the two-body slope parameter, and a is the ratio of the 
real part to the imaginary part of /nn(q — 0). For the inelastic terms, we write in the optical limit (ref. 13) 

Q] ( b/ ~ s j) Qj ( b - s j) = {exp [O (b' - s' , b - Sj )] - l} (32) 


with 


n 


(b' ~ s' , b - Sj) = ^ ' 2 E/ <iVVe" ) ' (b - ! i +s ">e ' q t b ' S i +S “)/NN(q)4 N (q') 


(33) 


Model for Nuclear Density Matrix 

We next describe a local density approximation for the one-body density matrix. For a projectile nucleus 
with Ap nucleons, the one-body density matrix is defined in terms of the complete nuclear wave function (4/) 
as (ref. 14) 

Kr,r') = y * 2 rfr 3 ,...,dr /lp 'I>t(r',r 2 ) ...,r j 4 fi )$(r,r 2 ,... > r / 4 p ) (34) 


The evaluation of equation (34) requires knowledge of the complete nuclear wave function; however, in practice 
a model is introduced. In the Fermi gas model (ref. 14), 


p(r,r') = p 0 


Sjlfejr-r 7 ;) 
kp\r — r'| 


(35) 


where k F is the Fermi momentum and p 0 is the density of nuclear matter. The Fermi gas model is known to 
provide a poor representation of the density matrix; however, its form suggests the use of a local density model 
where the density matrix factors are given as 


p(r, r') « p(y) n(x) (36) 

with x = r — r' and y = |(r + r 7 ), and where the one-body density is given by the diagonal part of the density 
matrix 


p(r) = p( r, r' = r) 

Here, n(x) is the Fourier transform of the nucleon momentum distribution 

(37) 


n(x) — J dp e* p x n(p) 

(38) 

where n(p) is defined by 

n (p) = J dr dr e zp x p{ r, r') 

(39) 

with normalization 

1? 

co 

1! 

I— 1 

(40) 


The one-body density is reasonably well known from elastic electron scattering. The nucleon momentum 
distribution at small to modest values of p is known from inclusive inelastic electron scattering. For large 
values of p, the backward production of protons suggests that large enhancements occur because of correlation 
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effects in contract to predictions of single-particle models and that the enhancements are largely independent 
of the nuclear mass (refs. 15-17). Before introducing phenomenological parameterizations of n(p), we note 
that by using equations (36) and (33), equation (19) can be written as 


dNi 

dk 


Jd 2 xd 2 pe ! < k -L Pj - ) xi n( v /p2 +fc 2) 


X 


d 2 yp(y ) 


3 D(b.x,y) 



(41) 


where kj^ is the transverse part of k and kp is the longitudinal part of kg. Equation (41) shows that the 
production spectrum in the longitudinal direction corresponds closely to the internal momentum distribution 
independent of the collision dynamics in the Glauber model, whereas the transverse spectrum is modified by 
multiple scattering. 

Haneishi and Fujita (ref. 17) have introduced 


3 

ra(p) e~ p2 / 2 Pi (42) 

i = l 

where C{ and p-i are constants listed in table 1 . The last term in equation (42) is expected to directly reflect 
the nuclear correlations. In equation (42), p x is related to the Fermi momentum by p 1 = sJYJlkp. Values for 
kp from experiments of Moinz et al. (ref. 18) are listed in table 2. For we use a value for kp corresponding 
to the matter radius of 1.69 fm. 

Table 1. Parameters for Momentum Distribution Table 2. Experimental Determination of Fermi 

Model of Equation (42) Momentum for Several Nuclei 


i 

Pi, MeV/c 

Ci 

1 


1 

2 

\/6/5k p 

0.03 

3 

500 

0.003 (0.0002)“ 


“Value found empirically. 


Nucleus 

k p , MeV/c 

12 C 

184 

40 Ar 

251 

208pb 

265 


Amado (ref. 19) has considered many body correlations by solving the Schrodinger equation directly for the 
case of delta-function potentials where, for large values of p, he finds 

n(p) “ itSFri (43) 

An analysis of backward proton scattering supports the value 7 = 90 MeV/c independent of the nuclear mass. 
We continue equation (43) to small values of p by using 


n(p) = n 0 


e~P 2 M + d 0 - 


7P 


sinh(qp) 

and consider several values for do below. We note that the normalization in equation (44) is given as 


(44) 


no = 


( 2 t rpf) 


2W/2 , d 0 7r 5 7 3 


-1 


(2 ^ 


(45) 





log /?(/•>) 


and for n(x) we find that 


n(x) = n 0 (27 rp 2 ) 3 / 2 e P?^ 2 / 2 + ^-^sech 2 ^^^jtanh^- 7 ^) (46) 

In figures 3, 4, and 5, we compare the models of equations (42) and (44) for 12 C, 40 Ar, and 208 Pb, respectively. 
For display, n(p) has been normalized to unity rather than to the condition given by equation (40). For a pure 



Figure 3. Logarithm of internal momentum distribution for Figure 4. Logarithm of internal momentum distribution for 

12 C versus momentum for several models. 40 Ar versus momentum for several models. 



Figure 5. Logarithm of internal momentum distribution for 208 Pb versus momentum for several models. 
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Gaussian model (d 0 = 0 in eq. (44)), the momentum distribution is vanishingly small above about 3 fm" 1 . 
The model of equation (42) is shown using C 3 - 0.003 and a value of C 3 = 0.0002 determined empirically, as 
discussed below. Including the correlation term with no Gaussian component (dashed line) is seen to show a 
large depletion of strength at small values of p in order to preserve the normalization condition. The use of 
equation (44) with 0.01 < d 0 < 0.1 will provide a reasonable momentum distribution at small values of p and 
account for some correlation strength at large values of p. 


Results and Discussion 

In figure 6 we show calculations for the momen- 
tum distribution of produced protons in 12 C on 12 C 
reactions at an energy of 1.028.4 GeV in the forward 
direction. The first three terms in the series of equa- 
tion (25) are shown corrected for the production of 
protons only with the final-state interactions (FSI) 
neglected. The invariant momentum distribution is 
shown which is found by multiplying equation (25) 
by the energy of the secondary proton. The calcu- 
lations are made by neglecting the correlation term 
(the third term) in equation (42). The leading-order 
term is seen to dominate, but with important contri- 
butions from the higher terms. The shape of the dis- 
tribution for the higher order terms is similar to the 
leading term except at high momentum values where 
they fall off more rapidly. In figure 7 we compare 
calculations which include the final-state interactions 



Figure 6. Calculations of invariant momentum distribution for 
proton production from abrasion in 12 C- 12 C collisions at 
1.0284 GeV. Shown are first three collision terms without 
final-state interactions. 


using the same internal momentum distribution as 
in figure 6. The absorptive effect of the FSI leads 
to a decrease of about a factor of 2 in the momen- 
tum distribution. This large decrease indicates the 
importance of cascade effects of the knockouts with 
the projectile fragments which must account for the 
decrease in figure 7. 

In figure 8 we compare calculations for proton 
production from abrasion for several targets with 
the experimental data from references 20 and 21 for 
inclusive proton production. The calculations shown 
by the solid line include the correlation contribution 
in equation (42). However, the strength has been 
adjusted to C3 = 0.0002 from the suggested value 
°f C3 = 0.003 from Haneishi and Fujita (ref. 17) 
because their value was found to be much too large 
in the high-momentum region. The dashed lines in 
figure 8 neglect the correlation term (do = 0). A 



Figure 7. Calculations of invariant momentum distribution for 
proton production from abrasion in 12 C- 12 C collisions at 
1.0284 GeV with and without final-state interactions. 
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Figure 8. Comparison of calculations with experiments of references 20 and 21 for proton production for 12 C collisions on several 
target nuclei at 1.028^4 GeV. Calculations include final-state interactions with and without correlation effects as described in 
text. 


treatment of phase space that is more accurate than 
the Glauber model employed here may allow for 
a larger correlation term; however, we note that 
most of the cited calculations on backwards proton 
production (refs. 15-17) neglect the absorptive ef- 
fects that would cause some renormalization of their 
results. The agreement between calculations and ex- 
periments in figure 8 is fairly good except in the low- 
momentum region where other processes will con- 
tribute to the inclusive proton production. The 
calculations shown include the production from the 
target knockouts, but these make only a small con- 
tribution at large momenta. The largest differences 
in the high- momentum region occur for the *11 tar- 


get where we note that the experimental results are 
obtained indirectly by using a CH 2 target. 

In figure 9 we compare calculations with exper- 
iments for proton production in 12 C- 12 C reactions 
at 2.062A GeV. The agreement between calculations 
is also good as in the lower energy data for the 12 C 
projectiles. In figure 10 we repeat the comparison 
made in figure 9 shown in the rest frame of the pro- 
jectile. Also shown for comparison are plots of the 
internal momentum distribution times the proton en- 
ergy normalized to the data. The comparison made 
in figure 10 illustrates the modification of the mo- 
mentum distribution by the collision dynamics. 
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Figure 9. Comparison of calculations with experiments of references 20 and 21 for proton production in 12 C- 12 C collisions at 
2.062 A GeV. 



-.6 -.4 -.2 0 .2 .4 .6 

p, GeV/c 


Figure 10. Comparison of calculations with experiments of references 20 and 21 for proton production in 12 C- 12 C collisions at 
2.062,4 GeV in projectile rest frame. Also shown is plot of E n(p) for models discussed in test. 


In figure 11 we show calculations of proton production in Ar-KCl collisions at 0 .8 A GeV at a scattering 
angle of 10° in comparison with the data of Nagamiya et al. (ref. 22). For calculations, the target is assumed 
to be Ar, and the agreement is again good in the high-momentum region. At low values of momentum, other 
contributions to the inclusive proton production, as discussed in the introduction of this report, are present 
and have not been estimated here. Production mechanisms other than those from abrasion will become more 
dominant at larger scattering angles. Finally, in figure 12 we show comparisons with the data of reference 23 1 
for proton production at 5° from 40 Ar projectiles interacting with Be and Cu targets. 
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p, GeV/c 

Figure 11. Comparison of calculations with experiments of reference 22 for proton production at 10° in 40 Ar- 40 Ar collisions at 
0.8A GeV. Data are for Ar-KCl and show calculations both with and without correlation term and protons from target. 


1 A vital comparison was also made with data found in an unpublished report by V. Perez-Mendez et al. (LBL-7278, Dep. of 
Physics, Univ. of California, 1978). 
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p , GeV/c 

Figure 12. Comparison of proton production on 40 Ar Cu and 40 Ar-Be collisions with production at 5°. Data are taken from 
reference 23 for 1.8 A GeV beams. 

Concluding Remarks 

The production of nucleons from heavy ion abrasion has been formulated by using a multiple-scattering 
model. The resulting formalism allows for the physics of nuclear abrasion, which is important in describing 
heavy ion fragmentation cross sections to be directly related to the inclusive yields of nucleons. Comparisons of 
theory with experimental data were made for inclusive proton production at small angles, and good agreement 
was found. The model comparisons indicate an important role for final-state interaction effects and the internal 
momentum distribution of nucleons in describing nuclear fragmentation. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
August 25, 1994 
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Appendix 

Kinematical Phase Space and Inclusive Momentum Distribution 

In this appendix we discuss the relationship of the multiple production terms of nucleons abraded from 
the projectile to the kinematical phase space and also consider an approximation using the Glauber multiple- 
scattering series in which energy conservation is ignored entirely. 

The scattering amplitude for the heavy ion collision is related to the cross section by the phase space of 
each particle that appears m the final state. We consider inclusive reactions in which a nucleon originating 
m the projectile is measured. For simplicity, the final target state is not considered and we use closure on 

these states with a single momentum vector denoted by p to represent these states. The cross section is then 
determined by 

(2tt) 4 n 

-f- d p F - e mw* - E f ) ^ - p /> ii>ii 2 (ai) 


da = 


X 


n=\j-l 


ere (3 is the relative projectile-target velocity, F* represents the prefragments, n is the number of nucleons 
knocked out of the projectile m the overlap region with the target, and i and / denote the initial and final 
states, respectively. The prefragment will decay through particle emission if sufficient energy is available. To 
include the phase space of decay products of F*, we write 


dpp* — dpp j j dp, 


r~ 0 


(A2) 


where r denotes the ions (if any) emitted in ablation. In considering nucleon production from ablation, we 
would study p r . We use the momentum-conserving delta function in equation (Al) to eliminate p from 
equation (Al). F 

We next consider using energy conservation in equation (Al). Working in the projectile rest frame, we 
transform dp x to dq, where q - Pt “ Px is the total momentum transfer in the collision, and we use the 

transfer" 0 We^en find* fUnCtl ° n ‘ n equation ( A1 ) to eliminate dq L , where q L is the longitudinal momentum 

r 1 T n. 

da = (27r) 4 ^2 d<l T K 

X 


n d v r 


Lr=0 


where the phase space factor is defined as 

1 -dE f 1 

K - ~ 

P dg L P 


En^i 

n=lj=l 


l T /! 


(A3) 


E x E f * 


E F * o PT ~ q L ) + E x ypPj cos Oj - q L J 

The momentum distribution of nucleons from abrasion is then 

f] FT P p Tl 

~ = (2rr) 4 E /*fr E / fit d Pj II d PrK\T fi ? 

X n= 1 j = 2 r=0 

which corresponds closely to equation (3) if we make the replacement 


(A4) 


(A5) 


K = J S(Ei — E f ) J dE F * 6(Ei - E f ) 


(A6) 

If energy conservation is ignored entirely, we would have K = 1. By using the participant-spectator arrangement 
oi the Glauber series discussed previously, the inclusive momentum distribution is shown to become 




(A7) 
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where 


and 



(A8) 

Id 

(A9) 


The result of equation (A8) indicates that the inclusive momentum distribution from abrasion essentially 
ollows the shape of the leading-order term in the approximations discussed because the absorptive factors in 
equation (A7) change slowly with increasing n for Ap ^>1. 
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